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Abstract 

It is known tliat including vector mesons stabilizes the size of a Skyrmion without 
the need for a Skyrme term. This paper provides the first results for static multi- 
Skyrmions in such a theory. The rational map ansatz is used to investigate multi- 
Skyrmions in a theory which includes the lo vector meson and has no Skyrme term. 
Bound states with baryon numbers two, three and four are found, which have axial, 
tetrahedral and cubic symmetries, respectively. The results reveal a qualitative sim- 
ilarity with the standard Skyrme model with a Skyrme term and no vector mesons, 
suggesting that some features are universal and do not depend on the details of the 
theory. Setting the pion decay constant and meson masses to their experimental values 
leaves only a single free parameter in the model. Fixing this parameter, by equating 
the energy of the baryon number four Skyrmion to the He^ mass, yields reasonable 
results for other baryon numbers. 
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1 Introduction 



Skyrmions are topological solitons that model baryons within a nonlinear theory of pions, 
arising as a low energy effective theory from QCD in the limit of a large number of colours 
[Ij. The standard Skyrme model [2] includes only the pion degrees of freedom and the 
Lagrangian requires the inclusion of a Skyrme term, which is quartic in derivatives. The 
role of the Skyrme term is to balance the sigma model contribution and provide a scale for 
the soliton, as required by Derrick's theorem [3J. The Skyrme term has some drawbacks, for 
example, it makes the theory non-renormalizable and the classical dynamical field equations 
have potential instabilities associated with the loss of hyperbolicity [1]. 

Over twenty years ago it was shown [S] that by generalizing the nonlinear pion theory to 
include vector mesons, the size of a Skyrmion is stabilized without the need for a Skyrme 
term. The original work included only the uj meson but later extensions added other vector 
mesons, for example p mesons O [71 [8] . These investigations produced promising results and 
revealed some improvements over the standard Skyrme model, although all these studies 
were limited to the sector with baryon number one. 

The single Skyrmion solution is spherically symmetric and therefore it can be constructed 
by solving only ordinary differential equations; though even these must be solved numeri- 
cally. However, multi-Skyrmions are not spherically symmetric and therefore highly nonlin- 
ear partial differential equations in three-dimensional space must be solved to study baryon 
numbers greater than one. The substantial difficulties that need to be overcome mean that 
even today there are still no results available on static multi-Skyrmions in theories including 
vector mesons; though it has been demonstrated that a product ansatz allows well-separated 
Skyrmions to be placed in an attractive channel [9], and recently a Skyrme crystal has 
been investigated [10]. The purpose of the present paper is to provide the first results on 
static multi-Skyrmions in theories including vector mesons. Details are presented for baryon 
numbers from one to four, but the methods described are also applicable to larger baryon 
numbers. 

Substantial progress in both numerical and analytic approaches to standard Skyrmions 
means that multi-Skyrmions are now fairly well understood in the Skyrme model (for a 
review see [H]). The approach taken in this paper is to apply some of these techniques, in 
particular the rational map ansatz [12], to study multi-Skyrmions in a theory without the 
Skyrme term but including the u vector meson. Briefiy, the rational map ansatz will be 
applied to provide an approximation to the Skyrme field and the u field will be computed 
using an expansion in terms of symmetry adapted spherical harmonics. 

Bound states with baryon numbers two, three and four are found, which have axial, 
tetrahedral and cubic symmetries, respectively. The results reveal a qualitative similarity 
with the standard Skyrme model, suggesting that some features are universal and do not 
depend on the details of the theory. Setting the pion decay constant and meson masses 
to their experimental values leaves only a single free parameter in the model. Fixing this 
parameter, by equating the energy of the baryon number four Skyrmion to the He^ mass, 
yields reasonable results for other baryon numbers. 

Finally, recent developments in AdS/QCD have led to renewed interest in baryons as 
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Skyrmions [13]. In particular, these studies point to the importance of the inclusion of vector 
mesons and provide a string theory motivation for old ideas of vector meson dominance. 

2 Skyrmions and the u meson 

The Skyrme field U takes values in SU (2) and satisfies the boundary condition that U ^ 1 
as |x| — i> oo. It is related to the triplet of pion fields tt through the formula 

U = a + iTz -T, (2.1) 

where r denotes the triplet of Pauli matrices and + tt • tt = 1. 

Topological solitons arise because there is a conserved topological current 

5^ = ^^e'^''"^Tr(9^f/ dJJ dpU U^), (2.2) 
247r^ 

and the associated integer topological charge B = J d^x is the soliton number and is 
identified with baryon number. 

The theory for the Skyrme field coupled to the uj vector meson is given by the Lagrangian 
density [S] 

io o 4 Z 

(2.3) 

where = 186 MeV is the pion decay constant, m.,^ = 138 MeV is the pion mass and 
rriuj = 782 MeV is the u mass. The constant (3 can be related to the uj ^ 37t decay rate, 
however, this is enhanced by the resonance u —>■ p + tc which is not included in the current 
theory. Therefore the experimental data only provides an upper bound on jS, which is found 
to be /3 < 25.4 |5j. 

It is convenient to remove various constants in the above Lagrangian density by rescaling 
the u field as u ujFt, and using energy and length units of F'^/m^ and 1/m^ respectively. 
The Lagrangian density then becomes 

1 1 1 

C = -Tr(9^f/9^f/t) +— Tr(f/-l)--(9^^,-9,^^)(9^^'^-9^^^) + -^^^^+(7^^i?^ (2.4) 

where M = m-^/m^ = 0.176 and g = [3mui/ Ft,. 

Note that the theory (12. 4p does not contain a Skyrme term and has only one free param- 
eter, (7, once the meson masses are fixed to their experimental values. 

Several values for g have been studied and the results are qualitatively similar in all 
cases. Results will be presented for the two values g = 96.7 and g = 34.7, which are both 
consistent with the above upper bound for f3. The first value, which is adopted from now 
on until further notice, is taken from [5], and is obtained by allowing both g and to be 
free parameters whose values are obtained by fitting the masses of the nucleon and delta 
resonance to the energies of a spinning B = 1 Skyrmion. This approach results in a value 
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of that is 124 MeV and therefore a httle less than the experimental result. The second 
value of g is calculated in a later section by setting to the experimental value and fitting 
the energy of the 5 = 4 Skyrmion to the He'^ mass. 

Only static fields are considered in this paper, therefore the spatial components of the 
topological current vanish = 0. As the topological current provides the source term for 
the field then its spatial components can be set to zero, uJi = 0. For ease of notation, in 
the following we drop the subscript on the temporal component and write uj = ujq. 

Taking into account the above comments, the static energy associated with the La- 
grangian density (12.41) is given by 

r ( 1 1 1 \ 

F = y ( — Tr (9,f/9,[/t) + -^Tr (1 - f/) - -d.ud^uj - -uo^ - guoB^ \ d?x. (2.5) 

Note the negative signs associated with the energy in the uj field, as it is a temporal compo- 
nent. 

Variation of the above energy with respect to uj yields the field equation 

didiU -uj = gB^, (2.6) 

which is a linear equation for u with a source term proportional to the topological charge 
density. The boundary condition is that u vanishes at spatial infinity. 

The problem at hand is therefore one of constrained energy minimization. The Skyrme 
field U is obtained as the field configuration, with a given topological charge B, that mini- 
mizes the energy (12. 5p . with u determined uniquely by U as the solution of equation (12. 6p . 
In the following section an approximation technique is described and implemented to solve 
this problem. 

Using equation (12. 6p . together with an integration by parts, the energy (12. 5p can also be 
written as 

E = \ ( -Tr(a.[/9.t/t) + — Tr(l -U)- -^9^B^\ d'x, (2.7) 

which will be useful later. Alternatively, the same procedure can be used to remove the 
last term in expression (12. 5p in exchange for changing the two minus signs to plus signs in 
front of the u dependent contributions. This form makes it transparent that the additional 
contribution to the energy from the inclusion of the u field is non-negative. 



3 Constructing multi-Skyrmions 

Extensive numerical computations of the full nonlinear field equations of the standard 
Skyrme model, with massless pions, have determined the minimal energy Skyrmions with 
baryon numbers up to i? = 22 ^14]. These numerical results can be reproduced with a sur- 
prising accuracy using an approximation employing rational maps between Riemann spheres. 



known as the rational map ansatz ^2J. This ansatz is briefly reviewed below. 
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In terms of the standard spherical polar coordinates r, 6, 0, introduce the Riemann sphere 
coordinate z = e**^ tan(^/2), given by stereographic projection of the two-sphere. Let R{z) 
be a degree B rational map between Riemann spheres, that is, R = p/q where p and q are 
polynomials in z such that max[deg(p), deg(g)] = B, and p and q have no common factors. 
Given a rational map R{z) the ansatz for the Skyrme field is 



f/(r, z) 



exp 



1 + \R\ 



-\R\^ 
2R 



2R 

\R\^-l 



where /(r) is a real profile function which satisfies the boundary conditions /(O) 
/(oo)=0. 

Applying the ansatz (13.11) to the baryon density (12.21) produces the expression 
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(3.2) 



from which it is simple to verify that the topological charge is indeed 5, as a consequence 
of the boundary conditions on /(r) and the fact that R{z) has degree B. 

For the particular choice R = z, the ansatz (13.11) reduces to the standard hedgehog ansatz 
for a spherically symmetric B = 1 Skyrmion. For B > 1 this ansatz does not provide any 
exact solutions, but for suitable choices of rational maps and profile functions, it provides 
excellent approximations to the true multi-Skyrmion solutions of the standard Skyrme model. 
In particular, the ansatz reproduces the correct symmetries of the true solutions, and the 
energy of the ansatz typically only overestimates the true energy of a multi-Skyrmion by the 
order of a percent. 

Substituting the ansatz (13. ip into the energy of the standard Skyrme model leads to only 
one contribution which is sensitive to the properties of the rational map beyond its degree. 
This contribution originates entirely from the Skyrme term and involves the coefficient 



X= — 
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(3.3) 



Thus, within the rational map ansatz, the problem of finding the minimal energy Skyrmion 
in the standard Skyrme model reduces to the simpler problem of calculating the rational map 
which minimizes the function X. Given the minimizing rational map, or more precisely the 
associated value of X, the profile function can then be determined by minimizing an energy 
functional for /(r). 

The X minimizing rational map for 5 = 1 is the spherically symmetric map R = z and 
for 5 = 2 it is the axially symmetric map R = z^, in agreement with the fact that the 
minimal energy B = 2 Skyrmion is axially symmetric [151 [HI [13 • 



For any B > 1 the map R 



is axially symmetric, but it is not the X minimizing 



map for B > 2. For B = 3 and 5 = 4 the X minimizing maps are the unique maps with 
tetrahedral and cubic symmetry, respectively, and are given by 



R 



iz 



VSiz^ 



R 



+ 2V3iz^ 



- 2V3iz^ 



(3.4) 
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These symmetries again agree with those of the numerically computed Skyrmions in the 
standard Skyrme model [IS]. 

In the vector meson model there is no Skyrme term and hence it appears that the co- 
efficient X is not relevant. However, notice that if the Laplacian term is neglected in the uj 
field equation (12.61) then oj is simply equal to a negative constant times the baryon density. 
Substituting this approximation into the energy (12.71) gives that the interaction term is a 
positive constant times the square of the baryon density. With the rational map expression 
for the baryon density (13.21) the angular part of this energy is precisely X, and therefore this 
term arises in the vector meson model as a leading order contribution in a derivative expan- 
sion. Alternatively, a formal solution for uj can be written in terms of the Green's function 
for the massive Klein-Gordon equation. If this representation is applied to the interaction 
term in the energy then it becomes a non-local expression involving two factors of the baryon 
density. The angular part of this is therefore a non-local version of X, which reproduces X 
exactly in the limit in which the Green's function is replaced by a delta function. 

The above arguments suggest that the X-minimizing maps are also the appropriate maps 
for the vector meson model, and indeed it will be shown that the maps (13. 4p yield low 
energy bound states, in contrast to the maps R = and R = z^, for example. This is strong 
evidence that in the vector meson model the minimal energy Skyrmions for B = 3 and B = 4 
have the same Platonic symmetries as in the standard Skyrme model. Additional evidence 
is provided by studying a (2+l)-dimensional analogue of this problem, where numerical 
solutions of the exact static field equations reveal an amazing similarity between solitons in 
the Baby Skyrme model and its vector meson version \19\. 




The approach applied below involves exploiting the symmetry of the Skyrme field to 
solve the field equation for u by using an expansion in terms of symmetry adapted spherical 
harmonics. This aspect is reviewed in the following. 

Given a subgroup G C S0{3) of spatial rotations, introduce the symmetric harmonics 
Ki[6,(j)), as a set of real orthonormal functions, each of which is a linear combination of 
spherical harmonics Yjm(6',0), that is. 



where the coefficients aim are chosen so that each function Ki is invariant under the group 
G. Strictly speaking, the functions Ki should carry an additional index, since there may 
be more than one symmetric harmonic with a given value of /. However, for the specific 
calculations and symmetries discussed below, only the values < / < 10 will be required, 
and there is a unique symmetric harmonic for each / in this range; therefore for simplicity 
the additional index will be suppressed. 

As an example, if G = 5*0(2), then the axially symmetric harmonics are simply Ki = Y/o- 
The angular part of the rational map generated baryon density (13.21) is given by 




Ki = aimXi 



(3.5) 



m=—l 




(3.6) 
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where the normahzation is such that the integral of b over the two-sphere is equal to B. 

For a G-symmetric rational map the angular baryon density b can be expanded in terms 
of symmetric harmonics as 

b = J2biKi, (3.7) 

I 

where 60 = B/ (2^/71) from the chosen normalization. Similarly, uj can also be expanded in 
terms of symmetric harmonics as 



UJ 



where hi{r) are real profile functions which depend only on the radial coordinate r. The 
boundary conditions are that hi{oo) = and hi{0) = for / > with /io(0) = 0. 

The baryon density expression (13. 2p together with the expansions (13.71) and (13. 8p reduces 
the u field equation (12. 6p to a set of profile function equations 

, + ^, 2^?fsinV 
+ -\^ + j = --^^^i- (3-9) 

Substituting the rational map ansatz (13. ip and the expansions (13. 7p and (13. 8p into the energy 
(12. 5p and performing the angular integration results in 

E = E^ + ^ r{r^f' + 2B sin^ / + 2M^r\l - cos /) + %f' sin^ / J] khi}dr, (3.10) 
2 Jo ^ 

where E^ denotes a contribution to the energy which has no explicit / dependence. 
The equation for / obtained from the variation of the energy (13.100 is 

r + -/' - 4 2/ - sin / + hh[ = 0. (3.11) 

For functions hi and / which satisfy equations (13. 9p and (13. lip the energy expression 
(12. 7p can be used to rewrite the energy as 

E=^ Hir^f + 2B sin^ / + 2Mh\l - cos /) + sin^ / ^ bihi}dr, (3.12) 

2 JO TT^ ^ 

which is similar to (13.100 except that there is no longer a term which is independent of / 
and the coefficient in the final term has been halved. 

Given a G-symmetric rational map R{z) the angular baryon density (13. 6p can be calcu- 
lated and the coefficients bi in the expansion (13.70 computed. The profile functions hi and 
/ can then be found by numerically solving the equations (13. 9p and (13. lip , using a heat 
flow algorithm. Finally, these profile functions are used to determine the energy from the 
expression (I3.12p . As mentioned earlier, it is found that truncating the expansions at / = 10 
is sufficient to produce results of the required accuracy. 
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Figure 1: Profile functions /(r) for B = 1,2, 3, 4. The curves shift to the right with increasing 
B. 



The simplest example is the spherical B = 1 Skyrmion with rational map R = z. There 
is only one spherically symmetric harmonic Kq = Yqq = l/{2^/^T) and hence only the profile 
functions / and need to be calculated. In this case the rational map ansatz is exact and 
the above procedure reproduces the earlier result [5j. For g = 96.7 the energy is found to 
be El = 44.04 and the associated profile function / is plotted as the solid curve in Figure [H 
The spherically symmetric field u = h^Ko is presented as the solid curve in Figure O 

Next consider the axially symmetric B = 2 map R = z^. As mentioned earlier, the 
axially symmetric harmonics are Ki = Yiq. However, the baryon density has an additional 
symmetry under rotations by 180° around an axis orthogonal to the S0{2) symmetry axis, 
which implies that only symmetric harmonics with even / are needed. A computation of the 
angular baryon density produces the expansion coefficients 

b = 0.5647^0 - 0M3K2 + 0.1077^4 - 0.0267^6 + O.OOGKs - O.OOliTio + . . . (3.13) 

A numerical solution of the profile function equations (for / < 10) results in the energy 
E2 = 87.44 < 88.08 = 2Ei, and hence this is a state which is bound against the break-up 
into two single Skyrmions. Note that the binding energy is quite small, being less than 1%. 
However, the rational map ansatz is exact for a single Skyrmion and is only an approximation 
for multi-Skyrmions, therefore any computation of binding energies which compares to single 
Skyrmions will be a slight underestimate. 

The profile function /(r) is displayed in Figure [T] (the curves shift to the right with 
increasing B). Although the u field now has a non-trivial angular dependence, it is useful 
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Figure 2: The spherical average < > as a function of the radius r, for 5 = 1,2, 3, 4. The 
curves shift up and to the right with increasing B. 



to plot the spherically averaged value 



which is displayed in Figure [2] (the curves shift up and to the right with increasing B) . This 
reveals that although the u field is again large and negative at the origin it takes its most 
negative value at a positive radius. This is to be expected because the source term for u is 
the baryon density and this has a toroidal distribution. An isosurface plot of u is displayed as 
the top right image in Figure [3|, and the toroidal distribution is evident. Although isosurface 
plots of uj resemble the associated baryon density surfaces there are some differences, for 
example, the baryon density vanishes at the origin for i? > 1 whereas, as seen above, the u 
field at the origin is substantial. 

Axially symmetric fields for larger baryon numbers can be studied by applying the same 
procedure as above to the rational map R = z^. However, it is found that all axially sym- 
metric fields of this type are unbound against the break-up into B single Skyrmions. For 
example, for 5 = 3 the energy is calculated to be Ef'"'^' = 138.24 > 132.12 = 3Ei. Of course, 
since the rational map approximation has been used to compute this energy it is expected 
that the true energy of the axial B = 3 field is slightly less than the above value, but the 
correction is unlikely to be large enough to yield a bound state. This is strong evidence 
that the B = 3 Skyrmion is not axially symmetric, in agreement with the usual Skyrme 
model. As mentioned earlier, with a Skyrme term and no vector mesons the minimal energy 
Skyrmion with B = 3 has tetrahedral symmetry and is described by the first rational map 
in (13.41) . In the following it is shown that in the u meson theory a tetrahedral bound state 




(3.14) 
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Figure 3: Isosurfaces of for i? = 1, 2, 3, 4 (to scale). The isosurface values are given by 
B = l,uj = -0.44; B = 2,iu = -0.40; B = 3, u = -0.30; B = A, u = -0.28. 



also exists with baryon number three, which suggests that it is the minimal energy B = 3 
Skyrmion. 

To study Platonic Skyrmions the symmetric harmonics Ki are required, where G is one 
of the Platonic symmetry groups. The construction of Platonic harmonics using invariant 
generating polynomials was introduced by Bethe and collaborators [20l[2T], and this is briefly 
reviewed below for the tetrahedral case. 

The ring of tetrahedrally invariant homogeneous polynomials is generated by three poly- 
nomials of degrees two, three and four. 



To obtain a tetrahedral harmonic Ki, the first step is to determine the most general degree 
/ homogeneous polynomial, ki, that can be constructed as a linear combination of products 
of the polynomials fl3.15p . The coefficients in the polynomial ki are then fixed (up to an 
irrelevant overall factor) by requiring that ki solves Laplace's equation. 

Note that for / = 2, 3, 5 there is only one possible contribution to ki, namely p2, Ps and 
P2P3 respectively. Of these three polynomials only p3 satisfies Laplace's equation, hence there 
are no tetrahedral harmonics with / = 2 or / = 5, and obviously none with / = 1. 

Given the tetrahedrally invariant polynomial ki then ^//r' is a function of only the angular 
coordinates 6, and, with a suitable normalization, is the required tetrahedral harmonic Ki, 
which can be decomposed into spherical harmonics Yim- In addition to the above construction, 
Platonic harmonics may also be determined using a projector technique [22], but in either 



_ 2 I 2 I 2 
P2 — X-^^ -r X2 ~r X3, 



Ps = X1X2X3, 



— 4 _i_ 4 I 4 



(3.15) 
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case the results are 



?■ - 1 5 I /l4 

Ko = Yoo, Ks = -^(^32 - ^3-2), = ^^j ^44 + 14-4 + ^ ^Y,o 



5 



i^6 = ^ |n4 + n-4-\/^no| , K-r = z^ lY,,-Yr_e + )Jf^{Yr2-Yj_2) | , 



1 / IQ^ / / 7 /22 

= — ( Yss + Ys-s + 2y ^(>84 + I8-4) + ^>80 1 , (3.16) 



K - Y +Y + ^ (Y +Y 1 .r^Y 



Substituting the first rational map in (13.41) into the angular baryon density (13.61) and 
expanding in terms of tetrahedral harmonics produces 

h = OM6K0 - 0.5857^3 - 0.0707^4 - 0.1277^6 + 0.026/^7 + 0.0027^8 + 0.022/^9 + 0.006i\:io + . . . 

(3.17) 

Using the above coefficients the associated profile functions can be computed (/ is presented 
in Figure [T] and the spherical average < u; > is plotted in Figure [2]) and the energy found to 
be E3 = 126.63 < 131.48 = E2 + E1, and therefore a bound state. A tetrahedrally symmetric 
u isosurface is displayed as the lower left image in Figure [31 

Turning attention to -B = 4, the required rational map is the second map in (13. 4p and has 
cubic symmetry. The cubically symmetric harmonics are, of course, a subset of those with 
tetrahedral symmetry. The tetrahedral polynomials p2 and in (I3.15P are also invariant 
under the cubic group, but ^3 is not. However, = pi is invariant under the cubic group, 
and indeed P2,Pa,P6 are the generating polynomials. This implies that the required cubic 
harmonics are the tetrahedral harmonics Ki (I3.16P with even /. 

The expansion coefficients of the B = A angular baryon density are given by 

b = l.USKo - OmiKi - OmOKe + 0.0767^8 + 0.009Kio + . . . (3.18) 

and the results of computing the profile functions are again displayed in Figure[T]and Figure[2J 
The energy is calculated to be £^4 = 159.67, and this confirms that this configuration is 
bound against the break-up into all possible lower charge clusters. A cubically symmetric uj 
isosurface is displayed as the lower right image in Figure [31 

The energy results for this coupling, which recall is g = 96.7, are summarized in the 
second and third columns of Table [H where the energy is also presented as a ratio to that 
of a single Skyrmion (this is convenient for the bound state comparison). 

A range of values for the coupling g have been investigated and the results are quali- 
tatively similar. As discussed earlier, the value g = 96.7 is taken from |5], which involves 
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Eb in MeV 


Experiment 


1 


44.04 


1.000 


22.53 


1.000 


996 


939 


2 


87.44 


1.985 


45.20 


2.006 


1999 


1876 


3 


126.63 


2.875 


65.88 


2.924 


2913 


2809 


4 


159.67 


3.626 


84.28 


3.741 


3727 


3727 



Table 1: The energy, Eb-, of the charge B Skyrmion and the ratio Eb/Ei, of the energy to 
that of a single Skyrmion, for two different values of the coupling g. For the second value of 
the coupling the energy is also given in MeV and the corresponding experimental value is 
listed for comparison. 

fitting properties of the single Skyrmion to the nucleon and delta, whilst treating F^^ as a free 
parameter. However, it is known that this method of parameter fitting can be problematic, 
in any Skyrme model, because of difficulties associated with the rigid rotor approximation 
[23] . Access to the properties of multi-Skyrmions allows the following alternative procedure 
to be applied. 

All the physical parameters, meson masses mTr, m^^, and the pion decay constant F^, 
are set to the experimental values listed earlier, leaving only the single parameter g to be 
determined. A sensible way to fit this parameter is to match the energy of the B = 4 
Skyrmion to the He^ mass, which is 3727 MeV. The reason this is a good approach is that 
the ground state of He^ has zero spin and isospin, therefore there are no issues to address 
regarding the inclusion of quantum energies associated with spin and isospin. Adopting 
this approach yields the value g = 34.7 and the associated Skyrmion energies, plus the 
experimental masses for comparison, are listed in columns four to seven of Table [H 

One point to note is that reducing g leads to a reduction in the relative binding energies. 
In particular, the data presented in Table [1] reveals that for this new value of the coupling 
E2 = 2Ei X 1.003, and therefore it appears that the B = 2 Skyrmion is not a bound state. 
However, recall that the rational map ansatz is only an approximation for B > 1 and the true 
energy is expected to be anything up to the order of a percent lower than the approximate 
value. As the value of E2 is so close to that of 2Ei then this correction is expected to 
produce a. B = 2 bound state, with a small binding energy. In fact, experimentally the 
deuteron binding energy is around 2 MeV, which is only about 0.1% of the deuteron mass. 
Therefore, if a reasonable comparison with experimental data is to be achieved, then even a 
tiny overestimate of the B = 2 energy would indeed be expected to result in an apparently 
unbound state. 

A comparison of the last two columns in Table [1] shows that all the calculated energies 
are within 7% of the experimental values, which is very reasonable given that the theory has 
only a single free parameter g. This accuracy is comparable to that in the standard Skyrme 
model, when the pion decay constant is treated as a free parameter whose value emerges as 
considerably lower than the experimental one. 

The Skyrmion energies presented are the classical energies and do not include quantum 
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contributions associated with spin and isospin. These quantum energies are traditionally 
added via a zero mode quantization involving the calculation of inertia tensors from the 
classical solution [21]. Expressions for the required inertia tensors have been calculated in 
terms of rational maps and profile functions [251 126] and therefore the results in this paper 
could be used to calculate these quantum corrections to the presented classical energies. 

4 Conclusion 

Multi-Skyrmions have been investigated in a theory without a Skyrme term but including 
the uj vector meson to provide a scale for the Skyrmion. The approach employed here 
involved using the rational map ansatz and therefore an obvious avenue for future research 
is to perform full field simulations of this theory, to test the accuracy of the approximations 
used. In particular, in the standard Skyrme model it is known that for baryon numbers above 
seven, there is an important qualitative difference between multi-Skyrmions in a theory with 
massive or massless pions [271 [2H] . It would be of interest to know if a similar result holds 
for the theory with vector mesons and no Skyrme term. 

Extending the current analysis to a theory including other vector mesons, starting with 
the p meson along the lines of [B], is clearly of importance. Furthermore, the emergence 
of a Skyrme model with vector mesons from AdS/QCD [13] suggests that a Skyrme term 
should also be included, even if it is not required to provide a scale for the Skyrmion. The 
results presented in this paper suggest that the vector meson and Skyrme term theories have 
very similar properties, and therefore one might expect that a theory which includes both 
contributions will also have similar universal features. 
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